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Abstract 

Background: Identification of protein-protein interactions is a fundannental aspect of understanding protein 
function. A commonly used method for identifying protein interactions is the yeast two-hybrid system. 

Results: Here we describe the application of next-generation sequencing to yeast two-hybrid interaction screens 
and develop Quantitative Interactor Screen Sequencing (QlS-Seq). QlS-Seq provides a quantitative measurement of 
enrichment for each interactor relative to its frequency in the library as well as its general stickiness (non-specific 
binding). The QlS-Seq approach is scalable and can be used with any yeast two-hybrid screen and with any next- 
generation sequencing platform. The quantitative nature of QlS-Seq data make it amenable to statistical evaluation, 
and importantly, facilitates the standardization of experimental design, data collection, and data analysis. We 
applied QlS-Seq to identify the Arabidopsis tholiono 1\/1L02 protein as a target of the Pseudomonas syringae type III 
secreted effector protein HopZ2. We validate the interaction between HopZ2 and 1\/1L02 in planta and show that 
the interaction is required for HopZ2-associated virulence. 

Conclusions: We demonstrate that QlS-Seq is a high-throughput quantitative interactor screen and validate 1\/1L02 
as an interactor and novel virulence target of the P. syringae type III secreted effector HopZ2. 

Keywords: Next-generation sequencing, yeast two-hybrid, high-throughput screening, Arabidopsis, Pseudonnonas 
syringae, type III effector, ML02, HopZ 



Background 

The Gram-negative bacterial pathogen Pseudomonas syr- 
ingae uses a type III secretion system (T3SS) to translo- 
cate type III effector (T3SE) proteins into the cytoplasm 
of plant cells. The primary function of these T3SEs is 
believed to be the suppression of plant immunity [1-5]. 
Some plant hosts are able to respond to this challenge 
via effector-triggered immunity (ETI), a defense 
response that is elicited when a plant resistance (R) pro- 
tein recognizes a specific effector protein either through 
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direct interaction, or indirectly via the action of the 
T3SE on its host targets [6,7]. The pathogen may 
respond by acquiring a new effector protein to suppress 
this recognition or by diversif)^ing away from recognition 
[7,8]. Thus, the pathogen and host each endeavor to 
gain the upper hand, resulting in a co-evolutionary arms 
race. 

There are -60 T3SE families identified in P, syringae, 
yet a majority of these remain functionally uncharacter- 
ized. A key to ascribing functions to these virulence pro- 
teins will be the identification of their host target 
proteins. In addition, since many T3SEs have evolved to 
suppress plant immunity, they can be used as probes to 
identify important components of resistance signaling 
pathways. 

Ltd. Tills is an Open Access article distributed under the terms of the Creative Commons 
s.org/licenses/by/2.0), which permits unrestricted use, distribution, and reproduction in 
properly cited. 
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The HopZ family of T3SE proteins is an evolutionary 
diverse family that is part of the YopJ T3SE superfamily 
found in animal and plant pathogens [9,10]. The HopZ 
family of P. syringae is composed of three distinct allele 
families (HopZl, HopZ2 and HopZ3), while HopZl also 
has three closely-related allele sub-families (HopZla, 
HopZlb and HopZlc). HopZla is most similar to the 
ancestral HopZ allele and is recognized by the ZARl 
resistance protein in Arabidopsis [9,11,12]. Although 
closely related to HopZla, HopZlb is only weakly 
recognized, and HopZlc is not recognized in Arabidop- 
sis [11]. HopZ2 and HopZS are more similar to YopJ 
superfamily members found in other phytopathogens 
and were likely acquired by P. syringae via horizontal 
gene transfer [9]. Both HopZ2 and HopZS have been 
demonstrated to enhance P, syringae growth on Arabi- 
dopsis [11,13]. Overall, the HopZ family displays 
remarkable functional diversification in Arabidopsis with 
members able to enhance bacterial virulence while 
others trigger ETI. Therefore, the targets of this T3SE 
family will likely include critical components of host 
immunity. 

The yeast two-hybrid (Y2H) system is a powerful tool 
to query protein-protein interactions [14,15]. Although 
several modifications of this method have been devel- 
oped, they all involve using a bait protein of interest to 
identify interacting prey proteins. In general, this can be 
done by using a bait to systematically test specific prey 
clones, or alternatively using a bait to identify interact- 
ing proteins from a pooled library of prey clones. The 
former method has been applied extensively in high- 
throughput fashion to generate high quality protein-pro- 
tein interactome maps [16,17]; nevertheless, the cover- 
age of these interactome maps is relatively low and is 
typically limited to model organisms, which have high 
quality libraries of cloned open reading frames (ORFs). 
If an ORE library is not available, a widely used alterna- 
tive is to screen a cDNA library; however, this approach 
carries the ascertainment and representation biases asso- 
ciated with all cDNA library methods. Additionally, 
screening of these biased libraries is typically limited by 
the throughput that candidate interactors can be 
sequenced. 

Recently, Vidal and colleagues have estabUshed a fra- 
mework to generate and assess high-throughput Y2H 
screens and established the first protocol for next-gen- 
eration sequencing of protein-protein interactomes 
[18,19]. Their Stitch-seq method employs PGR to conca- 
tenate the sequences of putatively interacting bait and 
prey proteins so that they comprise a single amplicon 
for downstream next-generation sequencing [19]. PGR 
stitching is done via common priming sites located 
downstream of both the bait and prey sequences. Their 
protocol was specifically designed for Y2H assays using 



Gateway sequences and clones, but can be generalized 
to other vectors and a wide variety of interaction assays. 
While unquestionably promising, there are some poten- 
tial limitations to this approach. The first relates to the 
size of the Stitch-seq concatenated amplicon, which is 
substantially larger than read lengths produced by cur- 
rent next-generation genome sequencers. As the authors 
note, in principle this obstacle can be overcome as read 
lengths improve or through paired-end sequencing, but 
the short reads generated by many next-gen sequencers 
may prove difficult to associate with a specific gene 
when working with random cDNA libraries generated 
from organisms with limited genomic resources. 
Another potential limitation arises from the long lengths 
of the stitched PGR products, which encompass the bait 
ORE, a linker, and the prey ORE, and the need for two 
rounds of PGR, which may result in PGR biases that 
influence the recovery of candidate interactors. 

Here we describe Quantitative Interaction Screen 
Sequencing (QlS-Seq), which couples split-ubiquitin 
yeast two-hybrid screening with Illumina next-genera- 
tion sequencing to rapidly identify interacting partners 
of a bait of interest. We employed this high-throughput 
and quantitative interactor screen to identify host targets 
of the HopZ family of T3SE proteins, and then demon- 
strate that these targets include components of plant 
innate immunity in plants. All members of the HopZ 
T3SE family (except for HopZ3) are membrane-asso- 
ciated [11,20], and as such, we used the split-ubiquitin 
yeast two-hybrid screen that utilizes a membrane-asso- 
ciated bait protein [21] in order to enrich for physiologi- 
cally relevant interactors. We used this approach to 
identify ML02 as an interactor of HopZ2, and con- 
firmed this interaction in vivo by bimolecular fluores- 
cence microscopy (BiEG). ML02 has a characterized 
role in powdery mildew resistance, but had not pre- 
viously been shown to contribute to P. syringae growth. 
We demonstrate that ML02 contributes to resistance 
against P, syringae in Arabidopsis and is required for 
HopZ2 virulence function. 

Results 

Evaluation of the cDNA prey library 

Our cDNA prey library was commercially made (Nor- 
clone Biotech Laboratories, Ontario) from RNA 
extracted from uninfected 4-5 week old Arabidopsis 
rosettes, as well as plants infiltrated with a virulent 
pathogen (P. syringae pv. tomato DG3000, P/:oDG3000), 
a non-virulent strain lacking the T3SS apparatus 
(PtoDG3000 AhrcG), and an avirulent strain translocat- 
ing a T3SE recognized by an Arabidopsis R protein 
(P^oDG3000 with AvrRpml, which is recognized by 
RPMl). Although it is common to amplify primary 
cDNA libraries after their initial construction, this step 
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can potentially introduce representation biases that may 
influence the interactor screen. We amplified the pri- 
mary library by semi-solid amplification as this method 
is believed to reduce overall amplification bias [22]. We 
first sequenced the primary and secondary cDNA prey 
libraries to assess representation and bias arising from 
the initial library construction and subsequent amplifica- 
tion. Amplified DNA was sequenced on an lUumina 
GA-IIx using standard protocols at the University of 
Toronto Centre for the Analysis of Genome Evolution 
and Function (CAGEF). Sequencing of the primary- and 
amplified libraries yielded 59.8 M and 5.8 M reads 
respectively, which read-mapped to -11 K Arabidopsis 
genes (4,119 M and 213 M bases of data. Additional file 
1). The cDNAs ranged from --40 nt up to --2,800 nt, 
with most not being full length since they were 



generated by random hexamer priming (Additional file 
2A, B). A scatterplot of the hits/locus for the two 
libraries revealed very high congruence (R^ = 0.96, Addi- 
tional file 2C), indicating that very little bias was intro- 
duced during the amplification process. 

To further evaluate the range of genes represented in 
the library, we analyzed the gene ontology (GO) annota- 
tions of the cDNAs recovered (Figure lA). Many biolo- 
gical processes were represented including metabolism, 
response to stress or stimulus, development, transport 
and signal transduction, although unknown was the 
most common (42%). The subcellular localization of 
many of the cDNAs was also unknown; however, loci 
were identified in virtually every cellular compartment 
including 13% associated with membranes. Genes with 
known molecular functions included those involved in 
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Figure 1 Characterization of cDNA library and HopZ2*"^''^ putative interactors. A. Percentage of primary cDNA library (left) and HopZ2^^^ 
interactors (right) encoding proteins belonging to gene ontology (GO) terms for biological processes, molecular functions and cellular 
components. * indicates categories that are missing for HopZ2^^^ interactors. B. Percentage of genes in library that are upregulated in response 
to biotic stress (from bacteria, oomycetes or elicitors of innate immunity). 
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protein binding, hydrolase activity, nucleic acid metabo- 
lism, transcription factors, transporters, and kinases. We 
further examined publicly available microarray data 
available through the CAGEF Bio-Array Resource (BAR, 
http://bar.utoronto.ca, [23]) to determine whether loci 
represented in our library were differentially regulated 
in responses to biotic stress from bacteria {P. syringae), 
oomycetes {Botrytis cinerea, Phytophthora infestans, 
Golonivomyces orontii), or elicitors of innate immunity 
(harpins, lipopolysaccharides, and an oomycete elicitor 
NPPl). 38% of genes were upregulated more than 2-fold 
in response to biotic stress while 46% did not respond 



to biotic stress (Figure IB). A further 16% did not have 
probe sets to detect transcriptional changes arising from 
biotic stress. None of our library loci were downregu- 
lated in response to biotic stress. 

Quantitative Interaction Screen Sequencing (QlS-Seq) 

Since most members of the HopZ family of T3SE pro- 
teins are membrane-associated by myristoylation [11], 
we adapted the split-ubiquitin yeast two-hybrid system 
that was developed for transmembrane bait proteins 
(Figure 2A and Additional file 3) [21]. Based on mem- 
brane-association studies of K-Ras, we constructed a 
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Enrichment = [(HopZ-Luciferase) / (HopZ+Luciferase+Library)]*100 



Figure 2 Schematic of QlS-Seq protocol. A. Outline of QlS-Seq screening protocol described in text. B. Example of enrichment calculation 
showing the effect of varying the recovery of a bait of interest (HopZ), the non-specific control (luciferase), and the frequency in the library. The 
counts and enrichment score are heat-mapped with high numbers in blue to low numbers in yellow. 
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bait vector with a C-terminal prenylation signal and 
polybasic sequence in order to stably associate the pro- 
tein with membranes (Additional file 4) [24,25]. 

We screened the Arabidopsis prey library by cloning 
catalytic mutants of all five of the HopZ alleles, as well 
as the Arabidopsis R protein ZARl that recognizes 
HopZla into our membrane-tethered bait vector (Figure 
2 A, Additional files 3, 4). We used the catalytic mutants 
of the HopZ T3SEs in order to prevent processing of 
putative substrates and potentially stabilize transient 
interactions that would occur with the active enzyme. 
This method has been termed inactive catalytic domain 
capture (ICDC) and has proven successful for stabilizing 
enzyme-substrate interactions [26,27]. We also screened 
the prey library with \io^Y2ptoi which is a P. syringae 
T3SE that is unrelated, both with respect to function 
and sequence, to the HopZ family, as well as luciferase 
as a non-specific protein control since it should not spe- 
cifically interact with any Arabidopsis proteins. We 
transformed our baits of interest by lithium acetate/PEG 
transformation and screened them on plates with drop- 
out media lacking histidine, one of the interaction 
reporter genes. Approximately 2000 colonies of each 
bait were collected and then restreaked twice on drop- 
out media lacking tryptophan in order to preferentially 
retain the prey plasmid and lose the bait plasmid 
thereby reducing the complexity of the plasmid pool for 
next-generation sequencing [28]. The colonies were har- 
vested en masse, digested with lyticase, and then plas- 
mids were purified using alkaline lysis. Finally, the prey- 
plasmid inserts were amplified with vector-specific pri- 
mers using low-cycle PGR to reduce amplification bias 
and the PGR product was Illumina sequenced. Each bait 
provided 4.7 M to 33.6 M quality reads (176 M to 2,544 
M bases of data. Additional File 1) which were read- 
mapped to the Arabidopsis reference using Novoalign 
(http://Novocraft.com), which performs base-quality 
aware, global alignments with affine gap penalties using 
full implementation of the Needleman-Wunsch algo- 
rithm. The number of reads per Arabidopsis coding 
sequence were converted to reads per million in order 
to normalize across samples. 

We then assessed for overall enrichment of each can- 
didate interactor by scaling the number of hits observed 
between our bait of interest and each candidate interac- 
tor relative to the frequency of those candidates in the 
primary library. The candidates' general 'stickiness' was 
also assessed by the number of times it was recovered 
using luciferase, our non-specific bait protein. Specifi- 
cally, enrichment was calculated as (T3SE-luciferase)/ 
(T3SE-i-luciferase+library) ''100, where each term is 
scaled as the number of mapped-reads per million (Fig- 
ure 2B). This enrichment measure scales from 0 to 100 
with higher values corresponding to those candidate 



interactors that do not bind luciferase (are not sticky) 
and are rare in the library (Additional file 5). 

Functional analysis of HopZ2-interactors 

We elected to focus our initial functional study on 
HopZ2 because it can promote P. syringae growth in 
Arabidopsis and also since a preliminary analysis of the 
data provided the most interesting candidate interactors. 
Our enrichment analysis of HopZ2^'^^^^ (hereafter 
HopZ2^^^) interactors identified several highly overre- 
presented and specific candidate interactors (Table 1 
and Additional file 5). HopZ2^^^ interactors were 
enriched for membrane-associated proteins (28% 
HopZ2*^^^ vs. 13% cDNA library; Figure lA) as well as 
proteins associated with responses to stress relative to 
the prey library (17% HopZ2*^^^ vs. 7% cDNA library; 
Figure lA). Based on our sequencing of the cDNA prey 
library, we could also assess the percent cDNA coverage 
of each HopZ2^^^ interactor in the prey library. For 
functional analyses we focused on the HopZ2^^^ inter- 
actors that had: (1) an enrichment value > 90% (33 loci), 
(2) were represented by clones > 33 amino acids in the 
cDNA prey library (ie. cDNA > 100 nucleotides) (18/33 
loci) and (3) were specific to HopZ2^^^ (i.e. not present 
with an enrichment score of > 50% in other baits tested) 
(11/33 loci). We hypothesized that these genes would 
include HopZ2 virulence targets and that their disrup- 
tion would alter P, syringae growth. 

We measured P, syringae growth in Arabidopsis lines 
carrying T-DNA insertions for each HopZ2^^'^ specific 
interactor to determine if the candidate HopZ2 interact- 
ing proteins played any role in P. syringae disease or 
resistance. We focused on interactors for which there 
were confirmed homozygous T-DNA insertion lines 
available and that were predicted to have the T-DNA 
insertion in an exon of the gene, and thus be more likely 
to interrupt the protein (5/33 loci. Table 1). We assayed 
for changes in immunity by infiltrating the T-DNA inser- 
tion lines with the virulent pathogen P^oDG3000 and 
evaluating bacterial growth after three days. Insertions in 
genes At4g35750, At5g20700, At4g00430 and Atlg68440, 
showed no difference in PtoDG3000 growth compared to 
the wild-type Gol-0; however, an interruption in gene 
Atlgll310 (line mlo2-7), encoding ML02 showed a ten- 
fold decrease in P^oDG3000 growth (Figure 3A). To 
further assess if this locus plays a role in resistance of 
Arabidopsis to P^oDG3000 we tested an additional T- 
DNA insertion Une in Atlgll310 {mlo2-6, Figure 3B). 
Bacterial growth was reduced by approximately 10-fold 
in mlo2-6 compared to Gol-0 wildtype (Additional file 
6A) indicating that the mlo2 mutation increases resis- 
tance to PtoDCSOOO. 

The mlo2-ll (pmr2-l) mutant of Arabidopsis, which is 
a point mutant (D287N) in ML02 Figure 3B), was 
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Table 1 Top HopZ2*^^'^ Interactors 
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^ rpm = reads per million 

^ Enrichment = [(HopZ-Luciferase)/(HopZ+Luciferase+Library)]^100 
^ See Additional file 5 for enrichment data with other baits. 

"Short" indicates that the prey cDNA library clone was less 100 nt in length. ND = not determined. 
^ "Exon" & "5' UTR" indicate that the T-DNA insertion is in either the exon or the 5' UTR. 
*Tested for increased or decreased bacterial growth 



previously identified in a genetic screen for mutants 
showing enhanced resistance to powdery mildew 
[29,30]. However, unlike mlo2-6 and mlo2-7, but consis- 
tent with the previous report [29], mlo2-ll {pmr2-l) 
does not exhibit any significant change in P, syringae 
growth compared to Col-0 (Additional file 6B). 

HopZ2 targets Arabidopsis N\102 

We further examined whether knocking out ML02 
would affect the virulence advantage normally conferred 
by HopZ2 when carried by the weak Arabidopsis 



pathogen P. syringae pv. cilantro 0788-9 (Pc/0788-9; 
[11]). We infiltrated Arabidopsis Col-0 and mlo2-7 
leaves with Pc/0788-9 carrying either an empty vector 
(EV) or HopZ2. As previously observed [11], HopZ2 
conferred a virulence advantage to Pc/0788-9 in Col-0, 
but the HopZ2-associated virulence advantage was com- 
promised in the mlo2-7 line (Figure 3C), indicating 
ML02 is required for HopZ2-associated virulence. 

To determine whether HopZ2 and ML02 interact in 
planta, we created fusions between both HopZ2^^^ and 
ML02 to the N- or C-terminus of YFP (nYFP or cYFP) 
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Figure 3 Disruption of ML02 compromises P. syringae 
virulence. A. The virulent pathogen PtoDCSOOO was pressure- 
infiltrated into the leaves of Arobidopsis Col-0 or T-DNA insertion 
lines in putative interactors of HopZ2^^^. An insertion in AtlgllSlO 
{mlo2-7) results in increased resistance to PtoDCSOOO. * indicates 
significant difference from Col-0 by Fisher's Protected Least 
Significant Difference (PLSD) test. Error bars indicate one standard 
deviation of the mean. B. Schematic showing the ML02 gene and 
protein structure with the point of insertion for the T-DNA insertion 
lines and the point mutant. The left side shows the gene with 
introns represented as lines and exons represented as boxes. ML02 
has two splice forms which differ at the 3' end of the gene. The 
right side shows the transmembrane structure of the ML02 protein. 
mlo2-6 and mlo2-7 are T-DNA insertion lines while mlo2-n is a 



point mutant. MLO 



' indicates the clone identified in the cDNA 



library screening. C. The weak pathogen Pc/0788-9 carrying the 
empty vector (Ev) or HopZ2 was pressure-infiltrated into the leaves 
of Arobidopsis Col-0 or mlo2-7. Statistical significance was 
determined as stated in part A. 



in a glucocorticoid-inducible conditional expression vec- 
tor. We used a partial clone of ML02 beginning at 
amino acid residue 281 of the full-length protein and 
containing the 4^*^, 5^^, 6^^ and 7^^ transmembrane 
domains as well as the C-terminal cytosolic tail of the 
protein (ML02^^'^^^), corresponding to the fragment of 
the clone in our cDNA prey library (Figure 3B). 

We infiltrated equivalent optical densities of Agrobac- 
terium carrying HopZ2::nYFP, HopZ2^^^::nYFP or 
HopZlc::nYFP with ML02'^^'^^^::cYFP, as well as the 
reciprocal combination. We used HopZlc as a negative 
control because it did not interact with ML02 in our 
yeast two-hybrid screening (Additional file 5, data not 
shown). Protein expression was induced by spraying the 
plants with dexamethasone post-infiltration. 72 and 96 
hours after dexamethasone application we observed 
bright fluorescence in leaf sections co-infiltrated with 
HopZ2::nYFP or HopZ2^^^:nYFP and ML02^^'^^^:: 
cYFP, as well as the reciprocal combination (Figure 4; 
Additional file 7A). No fluorescence was observed with 
HopZlc::nYFP and ML02'^^'^^^::cYFP (or the reciprocal 
combination) at these time points. The interaction 
between ML02'^^'^^^ and HopZ2 localized to the 
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Figure 4 HopZ2*^^^ but not HopZlc, interacts with ML02^^'^^° 
in planta by bimolecular fluorescence microscopy. Agrobocterium 
carrying HopZ2^/^::nYFP or ML02^^"2^°::cYFP (or the reciprocal 
combination) were mixed at equivalent optical densities and 
pressure-infiltrated into the leaves of N. benthomiono. Expression of 
the proteins was induced by 20 \xW\ dexamethasone. Sections of leaf 
tissue were imaged with a Leica SP5 confocal scanning microscope 
72-96 hours post-induction. The close-up images in the second 
column show the reticulate pattern reminiscent of the endoplasmic 
reticulum. The scale bar indicates 100 |jm. 
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periphery of the cell suggestive of the plasma membrane 
as well as reticulate network reminiscent of the endoplas- 
mic reticulum (ER; Figure 4). This localization pattern 
was also observed when an ML02^^'^^^::YFP fusion was 
transiently expressed in N. henthamiana (Additional file 
7B). 

Discussion 

We developed QlS-Seq, a quantitative, high-throughput 
yeast two-hybrid screening protocol combined with Illu- 
mina next-generation sequencing, to identify putative 
interacting proteins with the HopZ family of type III 
effector proteins. QlS-Seq provides many significant 
advances over traditional interactor screens: (1) it elimi- 
nates the need to individually sequence each interacting 
colony while at the same time vastly increasing the 
number of candidates interrogated; (2) the results are 
quantitative and therefore amenable to statistical analy- 
sis; (3) the method explicitly evaluates the enrichment 
of each interactor relative to both its presence in the 
prey library as well as its general (non-specific) sticki- 
ness; (4) sequencing of the prey library provides a hereto 
unprecedented ability to evaluate the cDNA library for 
complexity and completeness; (5) it is amenable to any 
type of yeast two-hybrid screen; (6) it is amenable to 
any type of next-generation sequencing; (7) it is comple- 
tely scalable and therefore applicable to experiments run 
in a very small, multiplex format, to very large auto- 
mated, high-throughput screens; and (8) the quantitative 
nature of the data also enhances the method's 'portabil- 
ity' among laboratories. 

A number of these points are worth elaborating. The 
ability to interrogate putative interactors relative to their 
presences in the prey library (points 3 and 4) is particu- 
larly critical when not working with well-established 
model systems. One of the great benefits of next-genera- 
tion sequencing is the ability to more easily study non- 
model systems. By definition, these systems have few 
established genomic resources, such as normalized 
cDNA libraries. The in silico normalization provided by 
QlS-Seq facilitates the use of any prey library, regardless 
of its means of preparation. For example, tissue, cell, 
age or stage-specific libraries could be rapidly con- 
structed and tested without the need for tedious and 
sample consuming normalization steps. 

Portability in the context mentioned in point 8, means 
that standards can be established for experimental 
design, data collection, and data analysis, which will 
allow experimental results to be comparable among 
laboratories. Examples of such portability standards 
included the MIAME (Minimum Information About a 
Microarray Experiment) [31] and MIGS (Minimum 
Information about a Genome Sequence) specifications 
[32]. Another benefit of these standards is that it 



encourages the development of data repositories and 
meta-analysis tools such as the Bio-Array Resource [23] 
for microarray data. 

A potential criticism of QlS-Seq is its cost-effective- 
ness, since the cost of next-generation sequencing is not 
trivial. Currently, it cost between US$1000-US$4000 for 
a single channel of lUumina next-generation data 
(depending on the specific protocol and platform). 
While we sequenced to quite high coverage in this 
proof-of-principle study, this depth is not generally 
required, and we found that 5 million reads were more 
than adequate. Since the current lUumina HiSeq2000 
platform currently produces over 100 million reads per 
lane, it should be possible to multiplex as many as 20 
samples per channel. Importantly, bar-codes can be 
directly incorporated onto the primers used to amplify 
the prey-plasmid inserts, thereby permitting the pooling 
of independent samples prior to lUumina sample prep. 
Early pooling of bar-coded samples means that only one 
sample prep is required for all pooled samples, and con- 
sequently, while the cost for a single channel of Illumina 
data may be US$2500, the cost per sample (if multiplex- 
ing 20 samples per channel) would only be US$125. 
This price is substantially less than the cost for Sanger 
sequencing 100 clones, and the cost will only continue 
to drop as the next-generation sequencing technology 
improves. 

Since we had previously shown that HopZ2 confers a 
virulence advantage to P. syringae in Arabidopsis, we 
therefore rationalized that we could use HopZ2 as a 
probe for the identification of innate immunity compo- 
nents. By conducting QlS-seq screens on all members of 
the HopZ family, we were able to identify proteins that 
interact specifically with HopZ2. These HopZ2 interac- 
tors were enriched for membrane-associated proteins as 
well as proteins from genes induced during stress 
responses, including Arabidopsis ML02. The lack of 
interaction between ML02 and the other HopZ family 
members suggests that the HopZ family has diversified 
to target different host proteins. 

ML02 has seven transmembrane domains with an 
extracellular N-terminus and an intracellular C-terminus 
and is localized to the plasma membrane [33]. HopZ2 is 
normally present at the plasma membrane and would be 
ideally localized to interact with ML02 [11]. Our analy- 
sis identified a partial clone of ML02 starting just prior 
to the fourth transmembrane domain and including the 
entire intracellular C-terminus (ML02'^^'^^^; Figure 3B). 
Using ML02'^^'^^^ in BiFC analyses, we demonstrated 
that HopZ2 and ML02 interact directly in planta. How- 
ever, the observed fluorescence complementation loca- 
lized to a reticulate structure reminiscent of the ER as 
well as the plasma membrane. This localization was also 
observed with ML02^^'^^^::YFP (Additional file 7B), 
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suggesting that ML02'^^'^^^ may be partially 
mislocalized. 

The MLO gene was first identified by map-based clon- 
ing in barley from mutants that were resistant to the 
powdery mildew fungal pathogen Blumeria (formerly 
Erysiphe) graminis f. sp. hordei {Bgh) [34]. However, 
m/o-based resistance in crop species has been employed 
by plant breeders for decades [35]. As in barley, Arabi- 
dopsis mlo2 confers increased resistance to a powdery 
mildew fungal pathogen, Golonivomyces (formerly Ery- 
siphe) orontii [30]. However, it has been reported (with 
data not shown) that P, syringae growth did not signifi- 
cantly differ in mlo2 compared to Col-0 [29,30]. Vogel 
and colleagues [29] evaluated symptom production in 
mlo2 (originally called pmr2) point mutants when infil- 
trated or sprayed with P^oDC3000, while Consonni and 
colleagues [30] tested bacterial growth and symptom 
production from PtoDC3000 in a T-DNA insertion line 
(SAIL_878_H12; mlo2-S) that is inserted towards the 
end of the 6^^^ exon (Figure 3B). Our growth assays with 
P^oDC3000 in the mlo2-ll {pmr2-l) point mutant con- 
firmed that it did not exhibit increased resistance to P. 
syringae (Additional file 6B). While mlo2-ll {pmr2-l), 
which has a D287N point mutation, in the intracellular 
loop between the third and fourth transmembrane 
domains [30], confers increased resistance to powdery 
mildew, it does not appear to be sufficient to confer 
increased resistance to P. syringae, suggesting that 
ML02 differentially contributes to immunity against 
these distinct pathogens. In contrast, unlike Consonni 
and colleagues we did observe a significant decrease in 
PtoDC3000 bacterial growth in two independent Salk T- 
DNA lines in the ML02 gene {mlo2-6 and mlo2-7), 
although we did not test their T-DNA insertion line 
{mlo2-S). Our results suggest that ML02 negatively con- 
tributes to resistance against P. syringae in Arabidopsis, 
and are consistent with the proposed role of ML02 as a 
negative regulator of defenses against oomycete 
pathogens. 

Previous work has shown that MLO is relocalized to a 
lipid raft-like domain in the plasma membrane upon 
pathogen attack [36]. ML02 has also been shown to 
negatively regulate PEN 1- dependent vesicular trafficking 
to regions of the plasma membrane associated with 
pathogen entry [7,37,38]. PENl is a syntaxin that has 
been associated with aberrant non-host resistance to the 
fungal barley pathogen Bgh, and is likely part of a 
SNARE complex involved in vesicular trafficking of 
defense components [37,39]. When PENl is recruited to 
sites of pathogen attack, it contributes to the rapid for- 
mation of papillae, an important component of the 
innate immunity [38,39]. Our data in conjunction with 
the prior studies suggest that pathogens may stabilize 
ML02 or cause its accumulation at the plasma 



membrane in order to suppress PENl-dependent secre- 
tion of defense components. 

There is precedence in the literature for P, syringae 
T3SEs to target negative regulators of plant immunity. 
The absence of RIN4 in rin4 rps2 plants compromises 
plant immunity whereas RIN4 overexpression enhances 
immunity [40]. Interestingly, at least four unrelated 
T3SEs have been demonstrated to target RIN4, poten- 
tially to enhance its role as a negative regulator of plant 
immunity [40-44]. Plant vesicular trafficking pathways 
are also targeted by P, syringae T3SEs. The T3SE 
HopMl induces the degradation of AtMIN7, an ARE 
guanine exchange factor (GEE) that is involved in vesi- 
cular trafficking [45]. Similarly, HopZ2 may stabilize 
ML02 in order to prevent the secretion of defense com- 
ponents to the regions of pathogen attack. 

Conclusions 

Overall we have demonstrated that QlS-Seq provides a 
powerful new approach to identify protein interactions 
using next-generation sequencing. We used this 
approach to identify Arabidopsis ML02 as a virulence 
target of the P, syringae T3SE HopZ2. Since HopZ2 (as 
well as other P. syringae T3SEs) is membrane localized 
we used the split-ubiquitin yeast two-hybrid system for 
interaction screening [11,21,46]. However, QlS-Seq is 
applicable to any sequencing-based yeast two-hybrid 
screening method. Eurthermore, this approach can be 
applied to both ORE as well as cDNA libraries. 
Although we sequenced the interactors of individual 
baits separately, the use of barcodes will allow the 
sequencing of pooled baits while maintaining the asso- 
ciations between interacting pairs. This approach will 
increase the number of baits that can be screened per 
experiment and decrease the cost of screening individual 
baits. In addition, as the costs of next-generation 
sequencing experiments continue to drop, QlS-Seq pro- 
mises to become a cost-effective alternative to tradi- 
tional yeast two-hybrid screening approaches. 

Methods 

Bacterial strains and routine culture conditions 

Escherichia coli and Agrobacterium tumefaciens were 
grown in Luria-Bertani broth, and Pseudomonas syringae 
was grown in King's broth (KB). Antibiotics were used 
at the following concentrations: for E, coli, 50 (ig/mL 
kanamycin, 100 (ig/mL ampicillin; for A. tumefaciens: 
100 (ig/mL kanamycin, 100 (ig/mL rifamcipin; and for P. 
syringae: 50 (ig/mL kanamycin, 50 (ig/mL rifamcipin 
and 50 (ig/mL cycloheximide. 

Plant growth conditions 

Arabidopsis plants were grown with 9 h of light (-130 
microeinsteins m'^ s'^) and 16 h of darkness at 22°C in 
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Sunshine #1 soil supplemented with 20:20:20 fertilizer. 
Nicotiana benthamiana plants were grown with 9 h of 
light (-130 microeinsteins m'^ s'^) and 16 h of darkness 
at 22°C in Sunshine #1 soil supplemented with 20:20:20 
fertilizer and osmocute. 

Cloning 

Pfu polymerase (Fermentas) was used for all cloning and 
all constructs were confirmed by sequencing. For the 
split-ubiquitin constructs, bait genes were amplified by 
PGR to contain an in-frame HA epitope, a polybasic 
region (K6 or K8) and a CAAX box, as well as appropri- 
ate unique restriction sites. The bait-HA-K6-CAAX 
genes were cloned into the pBT3-N vector (Dualsystems 
Biotech) using Sfil, The bait-HA-K8-CAAX genes were 
cloned into the pTLB-1 vector (gift of Dr. Igor Stagljar, 
University of Toronto) using Ncol, The orientation of 
each gene in the vector was confirmed. 

For the split- YFP constructs, the HopZ genes or the 3' 
end of ML02 were amplified by PGR to contain an in- 
frame HA epitope and appropriate unique restriction 
sites. All of the genes for the split-YFP system were 
cloned using Xhol and Stul into pBD-nYFP or pBD- 
cYFP. pBD-nYFP and pBD-cYFP were modified from 
pTA7002 [47] to add an HA tag and the N- or C-termi- 
nus of YFP between the Stul and Spel sites. The N-ter- 
minus of YFP encompasses residues 1-155 while the C- 
terminus of YFP includes residues 156 to the stop 
codon. 

The constructs used for plant infectivity assays were 
previously described [11]. In brief, the HopZ allele is 
expressed under its native promoter and contains an in- 
frame HA tag. 

cDNA library 

Five week old Arabidopsis rosette leaves were infiltrated 
by hand with a needleless syringe with P. syringae pv. 
tomato DC3000 (P/:oDC3000), P/:oDC3000 carrying 
AvrRpml or the AhrcC mutant of P^oDC3000 at an 
optical density of 0.1 (~5 x 10^ CFU/mL) at 600 nm. 
Infiltrated leaves were harvested at 4 hpi (P^oDC3000, 
P/:oDC3000 carrying AvrRpml, or P/:oDC3000 AhrcC) 
or 12 hpi (P/:oDC3000, P/:oDC3000 AhrcC). Uninfiltrated 
leaves were harvested at 4 pi and 12 hpi. RNA was 
extracted using Trizol (Invitrogen). mRNA was cloned 
into the pPR3-N vector (Dualsystems Biotech) using the 
Sfil sites (Norclone Biotech Laboratories, Ontario) with 
the NubG at the N-terminus of the prey proteins. The 
library contained -2.3 x 10^ clones, with an average size 
of 1.2 kB and was 90% recombinant. Amplification of 
the library was carried out by the semi-solid method 
[22]. 0.5 (iL of the primary library in E. coli strain DH5a 
was inoculated into 2x LB broth with 0.3% Seaprep 
agarose (FMC, Rockland) and 100 (ig/mL ampicillin. 



The inoculated cultures were then incubated in a water- 
ice bath for 1 hour. Subsequently, the inoculated cul- 
tures were incubated at 30°C for 44 hours without shak- 
ing. To sequence the primary and secondary libraries, 
low-cycle PCR amplification was carried out with pPR3- 
N vector-specific primers and a high-fidelity Taq/proof- 
reading polymerase mix (Fermentas, Burlington). This 
pool of DNA was sheared and prepared for lUumina 
sequencing by standard methods. 

Yeast two-hybrid screening 

HopZla^^^ HopZlb^^^ HopZlc^^^ HopZ2^^^ and 
ZARl^^ were expressed under the weak CYCl promo- 
ter in the pBT3-N vector, while HopZ3^^^ was 
expressed under the strong TEFl promoter in the 
pTLB-1 vector. AP-4 yeast [48] carrying the bait con- 
struct were transformed using the PEG/LiAc method. In 
brief, yeast carrying the bait construct were subcultured 
in 300 mL SD-Leu overnight to an optical density of 0.6 
at 600 nm. Yeast were washed twice in sterile H2O and 
resuspended in 1.5 mL. Transformations were per- 
formed with 1 (ig of cDNA library, 200 (iL of yeast cells 
and 600 (iL of PEG/LiAc (50% PEG, 120 mM LiAc, 10 
(iL 10 mg/mL boiled salmon sperm DNA) by the heat 
shock method at 42°C for 45 min. Yeast were washed 
twice in sterile H2O and plated on SD-LeuTrp and SD- 
LeuTrpHis + 3-amino-l,2,4-triazole (3-AT). Interacting 
colonies were identified by growth on SD-LeuTrpHis + 
3-AT. The appropriate amount of 3-AT was determined 
for each bait by testing for growth when transformed 
with the positive control pFur4-NubI and a lack of 
growth with pFur4-NubG [48]. Screening was performed 
until -2000 interacting colonies were identified. Colo- 
nies were restreaked twice on SD-Trp to preferentially 
lose the bait plasmid [10] and grown at 28°C. Prior to 
plasmid isolation, colonies were restreaked onto SD-Trp 
and grown at 28°C. Yeast were harvested en masse in 
SD-Trp and pelleted at 1000 g for 5 min. The pellet (- 
5 g) was washed in 0.1 M NaP04 pH 7.4 and 1.2 M sor- 
bitol and resuspended in 7.1 mL of lyticase buffer [0.1 
M NaP04 pH 7.4, 1.2 M sorbitol, 500 (iL lyticase 
(Sigma), 50 (iL of 10 mg/mL RNAseA]. 25 KU of lyti- 
case was resuspended in 0.01 M NaP04 pH 7.4 and 50% 
glycerol and used immediately. Yeast were incubated in 
the lyticase buffer overnight at 37°C. The Qiagen alka- 
line lysis spin kit was used to extract the plasmid DNA. 
The lyticase buffer was used instead of buffer PI. 
Volumes of buffers P2 and N3 were scaled up to the 
typical Qiagen proportions based on the volume of the 
resuspended pellet (1 vol lyticase buffer: 1 vol P2: 1.4 
vol N3). Lysis in the P2 buffer was done for 15 min at 
room temperature and 15 min at 65°C. Buffer N3 was 
chilled prior to use. After addition of buffer N3, the 
yeast were incubated on ice for 20 min. Yeast were 
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pelleted at 14000 rpm for 30 min at 4°C. The superna- 
tant was removed and cleared again by centrifugation at 
14000 rpm for 15 min at 4°C. The supernatant was 
loaded onto multiple Qiagen spin columns to purify the 
plasmid DNA. The columns were washed with buffer 
PB and buffer PE. Plasmid DNA was eluted with 50 (iL 
of buffer EB after a 1 min incubation. A second elution 
was performed with 35 \iL of buffer EB after a 1 min 
incubation. To sequence the putative interactors, low- 
cycle PGR amplification was carried out with pPR3-N 
vector-specific primers and a high-fidelity Taq/proof- 
reading polymerase mix (Fermentas, Burlington). This 
pool of DNA was sheared and prepared for lUumina 
sequencing by standard methods. 

Illumina Sequencing 

lUumina sequencing was performed either with 37 cycle 
single reads or 72 cycle paired-end reads (Additional file 
1) following the manufacturer s protocol on an Illumina 
GAIIx and pipelined using the GA pipeline vl.4. 

Bioinformatics 

Illumina reads were mapped to Arabidopsis gene models 
downloaded from NGBI, using a short read mapping 
tool novoalign (novocraft.com). From the mapping data, 
the number of mapped reads and the total length of 
mapped regions for each gene were determined with in 
house scripts. Read numbers per gene were further nor- 
malized as reads per million (rpm) within each sample 
and compared among the samples. The enrichment of a 
specific interactor with a bait of interest was determined 
by considering the number of reads obtained with the 
bait and luciferase and normalizing against the abun- 
dance of reads for luciferase, the bait and the library 
(Additional file 5). The percentage of the mapped length 
was calculated using length of mapped regions and the 
theoretical length of the gene model. Gene Ontology 
terms for Arabidopsis genes were downloaded from the 
TAIR website, and assigned to the genes in the cDNA 
library. Up- or down-regulation of each gene in 
response to biotic stress was determined from microar- 
ray data available through the CAGEF Bio-Array 
Resource (BAR, http://bar.utoronto.ca, [23]. Biotic stress 
treatments in the BAR included inoculation with viru- 
lent, avirulent and non-host P, syringae, inoculation 
with oomycetes {Botrytis cinerea, Phytophthora infestans, 
Golonivomyces orontii) and inoculation with elicitors of 
innate immunity (harpins, lipopolysaccharides, and an 
oomycete elicitor NPPl). 

Genotyping of T-DNA insertion lines 

The following T-DNA insertion lines were used: 
SALK_060284C (At5g20700), SALK_082464C 
(At4g35750), SALK_079850C (Atlgll310, mlo2-7), 



SALK_050191C (Atlgll310, mlo2-6), SALK_024490C 
(Atlg68440), SAIL_808_A10 (At4g00430). Genomic 
DNA was extracted from a leaf of 5-6 week old Arabi- 
dopsis plants. Primers were designed using the iSct fea- 
ture in the SIGnAL database. Primer sequences are 
available upon request. PCR-based genotyping was 
employed to determine the homozygosity or heterozyg- 
osity of the individuals. PGR products were sequenced 
using Big Dye Terminator 3.1 on an ABI 3730 genetic 
analyzer. 

Infectivity assays 

For infiltration, P, syringae was resuspended to an opti- 
cal density of 0.1 (~5 x 10^ GFU/mL) at 600 nm and 
diluted to a concentration of 1 x 10^ GFU/mL for 
growth curves. Plants were infiltrated by hand with a 
needleless syringe, as previously described [49]. Four 
disks (1 cm^) were harvested, ground in 10 mM MgGl2, 
and plated on KB with rifampicin and cycloheximide on 
days 0 and 3 for colony counting. 

Statistics 

For growth assays, 8-10 plants were used for day 3 
counts. Significance was determined using Fisher's Pro- 
tected Least Significant Difference (PLSD) on the day 3 
count data. Error bars indicate one standard deviation 
of the mean. 

Agrobacterium transient expression assays and BiFC 

Five-milliliter A. tumefaciens GV2260 cultures were 
grown overnight at 28°G in Luria-Bertani broth with 
kanamycin and rifampicin. The next day, the cultures 
were washed twice in induction medium (50 mM MES 
pH 5.6, 0.5% (w/v) glucose, 1.7 mM NaH2P04, 20 mM 
NH4GI, 1.2 mM MgS04, 2 mM KGl, 17 (iM FeS04, 70 
(iM GaGl2, 200 (iM acetosyringone) [50], and 3.75 mL 
was inoculated into 35 mL fresh induction medium to 
grow overnight. The following day, cultures were spun 
down, washed twice in 10 mM MES pH 5.6 with 200 
(iM acetosyringone and resuspended to an optical den- 
sity of 0.4 at 600 nm. The culture containing the 
ML02^^'^^^-cYFP plasmid was mixed in equal volumes 
with a culture containing the HopZlc-nYFP, HopZ2- 
nYFP or HopZ2^^^-nYFP plasmid. The culture contain- 
ing the ML02'^^'^^^-nYFP plasmid was mixed in equal 
volumes with a culture containing the HopZlc-cYFP, 
HopZ2-cYFP or HopZ2^^^-cYFP plasmid. The under- 
side of the leaves of 5- to 7-week-old N. benthamiana 
plants were infiltrated by hand with a needleless syringe. 
Plants were sprayed with 20 (iM dexamethasone (Sigma) 
1-2 hours after inoculation. Sections of leaves were 
imaged with a Leica SP5 microscope using Leica soft- 
ware at 24 hours (YEP fluorescence) or 72-96 hours 
(BiFG) post-dexamethasone induction. 
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Additional material 



Additional file 1: lllumina sequencing of cDNA libraries and 
interactors. Table showing the prey cDNA library used for the yeast 
two-hybrid screening, the number of lllumina cycles, the number of 
quality clusters, and the number of bases for each bait or the cDNA 
library. 

Additional file 2: Characterization of the Arabidopsis cDNA library 
used in this study. A. Histogram shows the number of Arabidopsis 
genes on log scale plotted against the length of the cDNAs present in 
the cDNA library. There are no genes present in the 2300-2499 or 2600- 
2699 ranges. There is one gene present in the 2500-2599 and 2700-2799 
ranges. B. Histogram shows the number of genes on log scale plotted 
against the percent coverage of the cDNAs in the cDNA library. C. 
Scatter plot of the genes present in the primary and amplified 
(secondary) libraries. The value of 0.9567 shows a high congruence 
between the primary unamplified and secondary amplified libraries. 

Additional file 3: Schematic of the split-ubiquitin yeast two hybrid 
approach. The N-terminus of the bait protein is compromised of the 
LexA-VP16 transcription factor (VP16), followed by the C-terminus of 
ubiquitin (Cub), the bait protein of interest, an HA epitope tag, a C- 
terminal polybasic region (K6 or K8) and CAAX box. We added the K6 (or 
K8) and CAAX box sequences to non-transmembrane bait proteins in 
order to anchor them to membranes and minimize autoactivation. At 
the N-terminus of the prey proteins is the N-terminal half of ubiquitin 
with an isoleucine to glycine mutation, which reduces non-specific 
association with Cub (NubG) [21]. The prey proteins may be membrane- 
associated or cytoplasmic (not shown). In the case of interaction 
between a bait and prey, the Cub and Nub portions of ubiquitin are 
brought into close proximity and reconstitute the full ubiquitin protein. A 
ubiquitin-specific protease (the "blades") cleaves the VP16 transcription 
factor, allowing it to translocate to the nucleus and activate the HIS 
reporter gene. 

Additional file 4: Bait and prey vectors used in QlS-Seq. The bait HA: 
K6:CAAX was cloned into the Sfi\ sites of the pBT3-N vector, in frame 
with an N-terminal LexA-VP16 (VP16) and the C-terminus of ubiquitin 
(Cub). Bait protein expression is driven by the weak CYCl promoter. 
Stronger bait expression was achieved by cloning the bait:HA:K8;CAAX 
into the Nco\ site of the pTLB-1 vector, in frame with an N-terminal VP16 
and Cub (not shown). The prey cDNAs were cloned into the Sfi\ sites of 
the pPR3-N vector, in frame with an N-terminal ubiquitin (the N-terminus 
of ubiquitin with an isoleucine to glycine mutation, NubG) [21]. Prey 
protein expression is driven by the weak CYCl promoter. 

Additional file 5: Rank list of top interactors with HopZ2 and their 
enrichment with other HopZ alleles, an unrelated T3SE HopF2 and 
ZAR1. Table indicates the percentage enrichment with each tested bait 
for the top HopZ2 interactors in Arabidopsis. 

Additional file 6: Disruption of IVIL02 compromises P. syringae 
virulence. A. The virulent pathogen ProDC3000 was pressure-infiltrated 
into the leaves of Arabidopsis Col-0 or the m!o2-6 T-DNA insertion line. * 
indicates significant difference from Col-0 by Fisher's PLSD test. Error bars 
indicate one standard deviation of the mean. B. The virulent pathogen 
ProDC3000 was pressure-infiltrated into the leaves of Arabidopsis Col-0 or 
the mlo2-n {pmr2-l) point mutant. Error bars indicate one standard 
deviation of the mean. 

Additional file 7: HopZ2 interacts with IVIL02^^"^^° in planta by 
bimolecular fluorescence microscopy and MLO^^"^^° localizes to a 
reticulate network. A. Agrobacterium carrying HopZ2::nYFP or ML02^^" 
^^°::cYFP were mixed at equivalent optical densities and pressure- 
infiltrated into the leaves of N. bentliamiana. Expression of the proteins 
was induced by 20 [iM dexamethasone. Sections of leaf tissue were 
imaged with a Leica SP5 confocal scanning microscope 72-96 hours 
post-induction. The scale bar indicates 100 |jm. B. Agrobacterium carrying 
I^^Q2Ai-28o..Ypp pressure-infiltrated into the leaves of N. 
benthamiana. Expression of the proteins was induced by 20 mM 
dexamethasone. Sections of leaf tissue were imaged with a Leica SP5 
confocal scanning microscope 24 hours post-induction. The scale bar 
indicates 100 pm. Fluorescence is observed in a reticular pattern 
reminiscent of the endoplasmic reticulum. 
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